Hyper novae /GRBs in the Galactic Center 
as possible sources of Galactic Positrons 
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ABSTRACT 

The observation of a strong and extended positron-electron line annihilation emission in the 
central regions of the Galaxy by INTEGRAL-SPl, consistent with the Galactic bulge geometry, 
without any counterpart in the gamma-ray range, neither at high energy nor in the 1809 keV Al 
decay line, is challenging. Leaving aside the geometrical question, we address the problem of the 
adequate positron sources, showing the potentiality of a new category of SN Ic, exemplified by 
SN2003dh, which is associated to a gamma-ray burst. This kind of supernova/hypernova/GRB 
event is interpreted as the result of a bipolar Wolf-Rayet explosion, which produces a large amount 
of 56 Ni and ejects it at high velocity along the rotation axis. The bulk of positrons resulting from 
56 Co decay escapes in the surrounding medium due to the rapid thinning of the cjccta in the polar 
direction. We show that a rate of about 0.02 SN2003dh-like events per century in the central 
region of the Galaxy is sufficient to explain the positron flux detected by INTEGRAL-SPl. In 
order to explain this flux by SN la events alone, a rate of 0.5 per century is necessary, much higher 
than indicated by Galactic evolutionary models applied to the bulge. Further observations of late 
light curves of SNe la and SNe Ic in the bulge of spiral galaxies, together with 3D hydrodynamic 
calculations of anisotropic ejections of 56 Ni in SN Ic/GRB events, will allow to estimate the 
separate contributions of SNe la and SNe Ic to positron injection. 



Subject headings: Galaxy: center 
gamma rays: theory 

Introduction 



gamma rays: bursts — supernovae: individual (SN2003dh) 



Since its discovery (Leventhal et al. 1978) the 
511 keV-line emission of the Galaxy is a prime 
target of gamma-ray spectroscopy (for reviews see 
Dermer and Murphy 2001, Milne et al. 2001, von 
Ballmoos et al. 2003). In the course of times, 
many instruments have been flown to detect this 
unambiguous signature of positron-electron anni- 
hilation, showing preferential emission in the gen- 
eral direction of the Galactic center. The mea- 
surement of the 511 keV line-profile by high- 
resolution germanium spectrometers has shown 
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that the line is essentially not shifted and narrow 
(about 2 keV FWHM). A step further in the un- 
derstanding of galactic positrons resulted from the 
CGRO-OSSE survey performed along the Galac- 
tic ridge, yielding a spectrum of the central regions 
of the Galaxy essentially consistent with that ex- 
pected from positronium annihilation (Kinzer et 
al. 2001). The OSSE observations have also sug- 
gested at least two emission components, a spher- 
ical one that can be linked to the Galactic bulge 
and a more elongated one that can be associated to 
the Galactic ridge (Milne et al. 2002). Indications 
of a third extended component situated at posi- 
tive Galactic latitudes have incited various spec- 
ulations about the underlying sources (von Ball- 
moos et al. 2003), in spite of the fact that its in- 
tensity and morphology is only poorly determined. 

Recently INTEGRAL-SPl has added its stone 
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to the construction, detecting a thin line (with a 
flux of ~f(T 3 ph cm" 2 s-\ Jean et al. 2003). 
This detection of a bright annihilation source in 
the central regions of the Galaxy confirms earlier 
findings on the line flux and width, but adds cru- 
cial information on the geometry of the source. 
Taking into account the spatial and temporal cod- 
ing by the instrument, Knodlseder et al. (2003) 
have provided new constraints on the morphology 
of the 511 keV e + e _ annihilation line emission 
from the central regions of the Galaxy. The emis- 
sion is concentrated in the direction of the Galactic 
center within a spherically symmetric zone whose 
FWHM extent is ^9 ° . The integrated measured 
intensity over such a feature leads to a flux of 
~10~ 3 ph cm" 2 s- 1 (Knodlseder et al. 2003). It 
is then firmly attested that the bulk of the 511 
keV annihilation emission seen in the central re- 
gion of the Galaxy originates in the Galactic bulge. 
Assuming that the bulk of the annihilation takes 
place at the distance of the Galactic center (8 kpc), 
via positronium, as implied by the OSSE observa- 
tions (Kinzer et al. 1996, 2001), the above re- 
ported flux of the 511 keV line emission infers a 
positron annihilation rate p = 7.7 x 10 42 (2-3//2) _1 
c + s _1 where / is the positronium fraction. Tak- 
ing / = 0.93 as derived from OSSE data (Kinzer 
et al. 2001), implies p = 1.3xl0 43 c+ s" 1 . In the 
steady state approximation (injection rate = an- 
nihilation rate), such an annihilation rate implies 
a similar positron injection rate. 

This high positron injection rate prompts us 
to reassess the question of the sources of the 
positrons. In the steady state approximation, only 
a restricted category of objects is, in principle, able 
to sustain such a high injection rate in the inner 
Galaxy, namely SNe la (Milne et al. 2001 and 
references therein), but just barely in the most 
optimistic case. This uncomfortable situation is 
a motivation to search for other potential sources 
of 56 Co, which, like SNe la, eject large amounts 
of 56 Ni, have a large expansion velocity, and are 
devoid of a hydrogen envelope that impedes e + 
escape. Such objects exist, they are the result of 
the explosion of bare C+O stellar cores (SNe Ic), 
i.e. exploding Wolf-Rayct (WR) stars. Unfortu- 
nately only a handful of them has been observed 
(Nomoto ct al. 2003, Maeda et al. 2003 and refer- 
ences therein) . These bright and short supernovae 
do not constitute a homogeneous class, and a case 



by case study is necessary. Here we focus on the 
most favorable example: SN2003dh, associated to 
the gamma-ray burst GRB030329 (Hjorth et al. 
2003, Stanek et al. 2003, Kawabata et al. 2003). 
In passing, we make a link between four aspects 
of the Galactic ecology: 

1. The bright 511 keV source sitting in the 
central regions of the Galaxy. 

2. The discovery of a special class of high lu- 
minosity supernovae, with fast ejecta, interpreted 
as the result of exploding WR, two of them be- 
ing associated to gamma-ray bursts, SN1998bw 
and SN2003dh (Nomoto et al. 2003, Woosley and 
Heger 2003, Maeda and Nomoto 2003a,b and ref- 
erences therein). 

3. The Galactic center region, spanning a few 
hundred parsecs, shows the presence of many mas- 
sive stars (Figer 2000), among which three massive 
star clusters: the Central cluster, the Arches and 
the Quintuplet clusters. Although small in size, 
the region contains ~10% of the molecular hydro- 
gen and young stars in the Galaxy. It is a major 
site of formation of massive stars, where O stars 
and WR stars abound. 

4. The existence of a bipolar Galactic wind in 
the center of our Galaxy, which could be induced 
by a past starburst (Bland-Hawthorn and Cohen 
2003). 

Beyond their capability to produce and inject 
significant amounts of 26 Al (Prantzos and Diehl 
1996), a WR explosion seems to lead to extraor- 
dinary phenomena, as, for instance a certain class 
of gamma-ray bursts associated to a bright super- 
nova of type Ic (SN2003dh). In this context, in- 
spired by the work of Woosley and Heger (2003) 
and Nomoto et al. (2003) we give a new role to the 
explosions of fast rotating WR stars (hypernovae) , 
that of generous positron injectors. 

2. In search of the positron source 

As detailed in the previous section, the positron 
injection rate related to the observed 511 keV line- 
flux is high, and the question of its origin is open. 
There is no lack of potential sources of Galactic 
positrons, but rather the contrary, and the best is 
to proceed by elimination. 

High energy and exotic processes, such as 
proton-proton interaction at high energy, neu- 
tralino annihilation (Bertone et al. 2002), decay 
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of Kaluza-Klein gravitons (Hannestad and Raffelt 
2003, Casse et al. 2003), are excluded since the 
positron production should be accompanied by a 
comparable flux of high energy gamma rays, not 
observed by CGi? O-EGRET (Hunter et al. 1997). 

One can also disregard the possibility that ac- 
creting black holes/microquasars (Mirabel et al. 
1992) could be able to account for the required 
positron injection within the Galactic bulge. In- 
deed, the GRANAT-SIGMA coded aperture tele- 
scope detected a transient line emission centered 
at 410 keV in the spectrum of the black hole can- 
didate IE 1740.7-2942 (the "great annihilator" , 
Bouchet et al. 1991), located 50' away from the 
Galactic center. If interpreted in terms of a red- 
shifted e + e~ annihilation line, the line flux of 
9.5xl0~ 3 ph cm~ 2 s _1 implies a positron anni- 
hilation rate of 5xl0 43 e + s" 1 (assuming a dis- 
tance of 8 kpc). However, it should be stressed 
that the Galactic bulge does not enclose enough 
compact sources of this kind to account for the re- 
quired rate given that they manifest only during 
very short periods. 

Isolated pulsars are also possible compact 
sources of positrons. But neutron stars young 
enough to be the site of the pulsar phenomenon 
are absent from the Galactic bulge since it is con- 
stituted by a very old stellar population. It is 
also tempting to invoke the production of mildly 
relativistic pair plasma outflow by a single gamma- 
ray burst event in the Galactic center within the 
past million years (Purcell et al. 1997). However, 
the annihilation line that may result from such an 
event (Furlanetto and Loeb 2002) is much broader 
than the observed one, and it remains to be shown 
how such a single event can so quickly populate 
homogeneously the whole Galactic bulge. 

The most promising source of Galactic positrons 
are the radioactive ejecta produced by the various 
nucleosynthesis sites. Indeed, through their [3 + 
decays, radioactive isotopes are capable of releas- 
ing enormous amounts of positrons, which annihi- 
late with electrons in the surrounding medium. 

Galactic novae can be excluded as dominant 
positron injectors, given their small positron pro- 
duction rate, and the absence of the associated 
22 Na line in the spectrum of the central regions of 
the Galaxy (Leising et al. 1998). 

The well-observed radioactive isotope 26 Al 



could produce significant amounts of positrons. 
However the flux of the 26 Al decay- line regis- 
tered in the central radian (3xl0 -4 ph cm~ 2 s _1 , 
Knodlscder 1999) and its thin latitude distribu- 
tion indicate clearly that this nucleus is not the 
main provider of the positrons that annihilate in 
the central regions of the Galaxy (see e.g. Purcell 
et al. 1997). 

The main source of positrons in this region, as 
well as in the whole Galaxy, is expected to be 
56 Co, the decay product of 56 Ni, which, in turn 
decays into 56 Fe, producing 19% of the time a 
positron. 56 Ni is produced explosively in super- 
nova events. SNe II yield rather modest amounts 
of 56 Ni (~ 0.1 MO against - 0.6 M© for SNe la, 
model DD23C, Milne, The and Leising 2001) and 
are covered by a thick hydrogen envelope that im- 
pedes the release of positrons in the surrounding 
medium. Thus SNe la remain the best positron 
fountains. They would surely produce positrons 
through the decay of 56 Co but the fraction of them 
escaping safely from the ejecta is difficult to ascer- 
tain. In favorable magnetic configurations, a few 
percent of them could leak out (Chan and Lin- 
genfelter 1993, Ruiz-Lapuente and Spruit 1998), 
but if the magnetic field is tangled, essentially no 
positrons escape. More empirically, based on the 
idea that positron escape should modify percepti- 
bly the late light-curve, Milne et al. (2002), have 
suggested that ~ 8xl0 52 positrons escape from a 
typical SN la. In order to produce the observed 
positron injection rate by SN la events alone, a 
mean SN la rate of 0.5 per century in the Galac- 
tic bulge is required. This rate is much higher 
than the predicted rate of 0.03 SNe la per century 
and 0.07 per century, as indicated by the models of 
Matteucci et al. (1999) and Nakasato and Nomoto 
(2003) respectively. 

3. Injection of positrons by exploding WR 
stars of the SN2003dh class 

Asymmetric explosions of WR stars producing 
a jet are under scrutiny (Nomoto et al. 2003 and 
references therein, Woosley and Heger 2003) . The 
work is still in its exploratory phase, but already 
important results emerge concerning the possibil- 
ity of a huge positron release in the framework of 
such events. Indeed the escape of positrons and 
its rate of increase depend on the distribution of 
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Ni in the velocity-mass space and on the opti- 
cal depth of the ejecta, which deserves a hydrody- 
namic study. However, in this exploratory work, 
semi-quantitative arguments would suffice. 

For the sake of simplicity, we assume that i) 
the radioactive nuclei are located deep inside the 
ejected envelope, ii) that the ejecta are in free 
homologous expansion, with an expansion veloc- 
ity v(r,t) oc r/t, where r and t are the radial and 
time coordinates, and Hi) that the propagation of 
positrons is purely absorptive. Under these as- 
sumptions the total positron optical depth is sim- 
ply r = npR where k is the absorption coefficient 
(opacity) of positrons, p the mean density, and R 
the radius of the ejecta. For an expanding ho- 
mogeneous sphere, pR scales as M 2 /E, where E 
is the kinetic energy imparted to matter by the 
explosion in units of 10 51 erg, and M the ejected 
mass in units of M© (Maeda et al. 2003). Now, if 
the ejection is bipolar, the matter being ejected in 
two opposite cones with opening half angle 9, the 
same relation holds providing that M represents 
only the mass of matter ejected within the cone 
and E the related energy. 

The transparency time of positrons (i.e. the 
time at which the optical thickness becomes ~ 1) 
in turn scales as ME -1 / 2 . The criterion which 
allows to select the best source candidate is a 
low ejecta mass over energy ratio and a high 
mass of 56 Ni produced. Compared to SNe la, 
SNe Ic are less efficient in the case of spherical 
explosions since their ejected mass is in general 
higher than 1.4 M© (the mass of typical SN la 
ejecta). But the situation is completely different 
in the case of asymmetric WR explosions. Ex- 
plosions with marked asymmetries lead i) to high 
velocities along an axis. The higher velocity in 
the favored direction allows an earlier escape of 
photons (at all energies) and positrons, due to 
a faster thinning out of the material; ii) to a 
strong mixing of radioactive nuclei ( 56 Co). In- 
deed, such conditions are fulfilled by the models 
of Maeda and Nomoto (2003a) and Woosley and 
Heger (2003), hereafter WH. The efficiency of the 
energy deposition of gamma-rays and positrons 
would be lowered resulting in a quick decline of the 
light curve, as observed in the case of SN2003dh 
(WH). We take advantage of this effect to suggest 
that SN2003dh-like events could inject a copious 
amount of positrons in their surroundings. 



To explain the peculiar light curve of SN2003dh, 
especially short and intense, WH have made a two 
components model (slowly equatorial ejection plus 
fast polar ejection). The exploding object is a WR 
star, having lost its H and He envelope, and thus 
reduced to a pure C+O core whose mass is about 
10 MQ. Considering that a compact central rem- 
nant of about 2 M© is formed, the ejected mass 
within a cone with opening half angle of 45 ° is 
1.2 MO. The energy released in each cone, 1.25 
10 52 erg, has been derived from Fig. 6 of WH. 
In such a case, the ME -1 / 2 ratio is 0.34 com- 
pared to 1.4 for type la supcrnovae, leading for 
SN2003dh to a transparency time about 4 times 
shorter. From the peak of the light curve, the 
total amount of 56 Ni (0.5 M©) is consistent with 
that of SN1998bw, of similar brightness. 

Exploiting the fact that both SNe la and hy- 
pernovae concern the explosion of C+O cores, we 
calibrate the escape fraction of SN2003dh in par- 
ticular on that of the sample of SNe la studied 
by Milne, The and Leising (2001). Positron trans- 
port depends sensitively on the strength and con- 
figuration of the magnetic field, which is unknown 
(Chan and Lingcnfclter 1993, Ruiz-Lapuente and 
Spruit 1998). Three extreme cases have been sug- 
gested: i) weak field with no confinement, ii) very 
confining strong and tangled field, and Hi) strong 
field with radial field lines. Following Milne, The 
and Leising (2001) we have adopted the radial 
field geometry favored by the analysis of late SN 
la light curves. Comparing the observed V band 
light curves from 22 SNe la to a model-generated 
one, Milne, The and Leising (2001) conclude that 
the late light curve is best fit by radial escape of 
positrons after ~ 390 days, leading to a number 
of ~ 8xl0 52 positrons escaping (i.e. 3.3% of the 
total number of positrons produced). 

Scaling this transparency time with ME -1 / 2 
and taking a mass of 1.2 M© (O + C + Ni) ejected 
in a cone with an opening half angle of 45 ° (as WH 
assume), with a kinetic energy of 1.25xl0 52 erg, 
transparency to positrons is achieved early (about 
107 days), and thus a significant fraction (about 42 
%) of positrons produced by the decay of 0.5 M© 
of 56 Co are expected to be released in the environ- 
ment. This fraction is a lower limit, since dynam- 
ical instabilities in the fast moving bipolar ejecta 
(WH model) would bring the radioactive nuclei 
closer to the surface, as indicated by the short rise 
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time of the light curve. Each SN2003dh-like event 
would then release at least about 10 54 positrons 
and most probably up to 2xl0 54 in the case of 
efficient mixing. We are inclined to consider that 
this kind of event offers a providential source of 
w 1 MeV positrons. Of course more quantitative 
work is necessary to strengthen this conclusion. 
We urge specialists of stellar models to calculate 
the positron fraction leaking out from asymmet- 
ric SN Ic explosions (in general and SN2003dh in 
particular), with detailed hydro-codes. 

It is worth noting that Mazzali et al. (2003) in- 
terpret the light curve of SN2003dh in terms of a 
composite symmetrical model. Here again, a large 
amount of 56 Ni resides in high velocity ejecta. It 
would be worth pursuing the study of positron es- 
cape in this framework. Other hypernovae/SNe 
Ic showing wider light curves (as e.g. SN1997ef, 
SN1998bw and SN2002ap, Maeda et al. 2003) 
seem less favorable. 

Now, coming back to the WH model, if the oc- 
currence rate in the bulge is one SN2003dh hyper- 
nova/GRB for ~ 25 SNe la, the two sources have 
equal impact. This hypernova rate seems at first 
glance rather high, due to the exceptional nature 
of SN2003dh. However selection effects could be at 
work, both geometrical and photometrical. First, 
if the model is right, SN2003dh is a bipolar su- 
pernova seen pole on; bipolar SNe observed in the 
equatorial direction are expected to be fainter than 
SN2003dh (Hoflich, Wheeler and Wang 1999) and 
more likely escape detection. Second, SN2003dh 
is exceptionally bright compared to other core col- 
lapse supernovae, but not to a SN la. Thus it 
would be quite difficult to detect this kind of ob- 
jects against the bright background of the central 
regions of spiral galaxies (Barbon et al. 1999). 
In this vein, it is not surprising that SN2003dh is 
hosted by a star-burst irregular galaxy (Matheson 
et al. 2003). Third, if the Galactic center has 
experienced, a few million years ago, a burst of 
massive star formation, as indicated by the Galac- 
tic wind which seems to escape from it (Bland- 
Hawthorn and Cohen, 2003), the SN Ib,c rate has 
been exceptionally high in the Galactic center. 

Indeed a rigorous and consistent calculation of 
positron escape from all SN I types (la, lb, Ic) tak- 
ing into account internal mixing and anisotropic 
ejection of 56 Ni has to be performed before con- 
cluding. Even SNe la show some signs of asym- 



metry through their polarized light (Wang et al. 
2003). Low-mass SNe lb (like SN1994I), which 
have not been treated specifically (except in the 
work of Ruiz-Lapucntc and Spruit 1998), are es- 
pecially interesting. 

4. Conclusions 

We have presented semi-quantitative argu- 
ments in favor of a certain type of SNe Ic, reserv- 
ing to a forthcoming article a more detailed study 
(Lehoucq et al. in preparation). We only have 
attacked the problem of the number of positrons 
injected in the region, postponing the description 
of their propagation-annihilation in conditions ap- 
propriate to the Galactic bulge to a subsequent 
work (Paul et al. in preparation). We have ad- 
vocated that SN2003dh-like events, arising from 
asymmetric explosions of Wolf-Rayct stars, can 
compete with SNe la as leading positron sources. 

The number of positrons injected by one hy- 
pernova/GRB (2xl0 54 ) is ~ 25 times higher than 
that injected by a SN la (8xl0 52 ) and consider- 
ably higher than that produced by a microquasar. 
In order to explain the detected positron produc- 
tion rate in the Galactic center bulge (1.3xl0 43 
c+ s _1 ) solely by the explosion of SN2003dh-type 
hypernovae, their average occurrence must be of 
the order of 0.2 per millcnium. In the case of SNe 
la solely, the occurrence must be 0.5 per century, 
which is above the prediction of current models 
(Matteucci et al. 1999, Nakasato and Nomoto 
2003). 

INTEGRAL continues to collect information on 
the inner Galaxy. More data on the 511 keV emis- 
sion of the central regions will become available in 
the near future, including a more detailed map of 
the bulge and lower halo. If the emission proves 
to be patchy, discrete events of the kind proposed 
will be favored. On the contrary, if it is symmet- 
ric, and showing a smooth gradient outwards, a 
new form of dark matter (scalar and light, Boehm 
et al. 2003) could be promoted. 

Further observations of late light-curves of SNe 
la and SNe Ic in the bulge of spiral galaxies and 
starburst galaxies will allow to estimate the sepa- 
rate contribution of SNe la and SNe Ic to positron 
injection. The study of the late light-curves of SNe 
la and SNe Ic could be a by-product of the cos- 
mological supernova surveys performed on ground 
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(e.g. MEGACAM, Boulade et al. 2003) and fore- 
seen is space (SNAP, Lampton et al. 2002). 

We warmly thank the referee, Keiichi Maeda, 
for pertinent and constructive comments. We are 
profoundly indebted to Francesca Matteucci, Nao- 
hito Nakamura and Ken Nomoto for having made 
available their estimate of the SN la rate in the 
Galactic bulge. The authors would like to thank 
the team of SPI, the Spectrometer aboard the ESA 
gamma-ray space telescope INTEGRAL, and Lau- 
rent Vigroux for useful discussions. 
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